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Abstract: Echinococcus granulosus is the causative agent of cystic echinococcosis with medical and veterinary impor- 
tance in China. Our main objective was to discuss the genotypes and genetic diversity of £ granulosus present In domes- 
tic animals and humans In western China. A total of 45 hydatid cyst samples were collected from sheep, humans, and a 
yak and subjected to an analysis of the sequences of mitochondrial cytochrome b {cytb) gene. The amplified PCR product 
for all samples was a 1 ,068 bp band. The phylogenetic analysis showed that all 45 samples were identified as £ granulo- 
sus (genotype G1). Ten haplotypes were detected among the samples, with the main haplotype being HI. The haplotype 
diversity was 0.626, while the nucleotide diversity was 0.001. These results suggested that genetic diversity was low 
among our samples collected from the west of China based on cytb gene analysis. These findings may provide more infor- 
mation on molecular characteristics of £ granulosus from this Chinese region. 
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Cystic echinococcosis (CE, also known as hydatid disease), a 
globally important and pathogenic zoonoses, is caused by the 
atiological agent, Echinococcus granulosus. China is one of the 
most important endemic regions of CE [1]. £. granulosus is now 
considered as a complex consisting of at least 4 species includ- 
ing E. granulosus sensu stricto (G1-G3 genotypes], E. equinus 
(G4), and E. ortleppi (C5). There are controversies for consider- 
ing G6 to GIG as 1 or 2 species, namely E. canadensis and 
E.intermedius [2-4]. Among them, E. granulosus s.s. is known to 
have a broad geographical distribution and a wide host range 
[5]. G1-G3 genotypes are also named the sheep (Gl), theTas- 
manian sheep (G2), and the buffalo (G3) strains, respectively. 
In China, the sheep strain is the predominant epidemic strain, 
while G3 and G6 strains also exist [6-14]. 

Increased knowledge of the genetic stmcture of a parasite 
could provide more information about its epidemiology and 
transmission, since genetic variants may affect infectivity and 
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pathogenicity [15]. In China, some mitochondrial and nuclear 
genes have been applied in the analysis of genetic information 
in E. granulosus, including cytochrome c oxidase subunit 1 
[coxl), NADH dehydrogenase 1 [nadl], ATP synthase FO sub- 
unit 6 [atp6), cytochrome h [cyth), elongation factor-1 alpha 
[efla), 12S, and 16S [2,13,14,16,17]. However, the genetic di- 
versity of E. gmnulosus in China has never been characterized 
based on the complete cytb gene, and only the partial cytb gene 
have been reported [17,18]. Previous studies demonstrated 
that both the complete cytb and coxl sequences showed 1 
causative agent of CE in a cat from Russsia, i.e., E. granulosus s.s. 
(Gl) [19]. Therefore, it is feasible that we use the complete qtb 
gene as a genetic marker to analyse genetic information of E. 
granulosus. The main purpose of this study was to figure out 
the genotypes and genetic diversity of E. granulosus present in 
westem China. 

A total of 45 hydatid cysts were obtained fi-om several hosts 
(sheep, yak, and humans) in Tibet Autonomous Region, Qing- 
hai, and Sichuan Province, China. Table 1 summarizes the 
number and origin of isolates examined in the 3 localities. All 
samples were stored at -70°C. Genomic DNA was extracted 
from each cyst by proteinase K treatment and phenol-chloro- 
form extraction. All DNA samples were stored at -20°C. 



205 



206 Korean J Parasitol Vol. 52, No. 2: 205-209, April 201 4 



Table 1 . Geographic origin and liost of China isolates of Echino- 
coccus granulosus identified by DNA sequencing of mitochondri- 
al cytb gene 



Province 




Host 




Total 


Sheep 


Human 


Yak 


Qinghai 


28 


0 


0 


28 


Sichuan 


0 


9 


1 


10 


Tibet 


7 


0 


0 


7 


Total 


35 


9 


1 


45 




Based on the published complete mitochondrial genome of 
E. granulosus (GenBank accession no. AF297617), the cytb gene 
was amplified using the following primer pairs: forward 
5'-AGATATrAGAGATGGAGGAG-3' and reverse 5'-AGCGGA 
GITAATACAGTCAG-3'. The primers were designed by Primer 
5.0 software, and the produa was 1,609 bp that contained the 
complete sequence of cytb gene. PGR was performed in a 25 pi 
reaction mixture containing 1 pi of template (the genomic DNA 
of each sample was used as a template for PGR), 1 pi of each 
primer (10 pmol/pl), 12.5 pi of 2 x Taq PGR MasterMix and 9.5 
pi of ddHiO. Negative controls were always included in PGR re- 
actions to assess possible contaminarion. The PGR protocol 
consisted of 30 thermal cycles with a denaturation step (90 sec 
at 94°G), a hybridization step (90 sec at 46°G) and an elonga- 
tion step (90 sec at 72°G) for each cycle. Afterwards, 5 pi of the 
PGR product was run in a 1% agarose gel and stained with 
ethidium bromide to detect the PGR amplicons. Amplification 
products were purified by using the Universal DNA Purification 
Kit (Cowin Biotech, Beijing Ghina) and subjected to automat- 
ed sequencing following subcloning into the pMD19-T vector 
(TaKaRa, Beijing, Ghina). All amplicons were sequenced at least 
3 times. 

Multiple alignments of nucleotide sequences were performed 
using the Clustal X 1.83 software, with variation sites deleted 
fi-om the alignments. Population diversity indices (nucleotide 
diversity number of haplotypes, and haplotype diversity) were 
calculated using DnaSP 5.0 software [20]. Phylogenetic trees 
were constructed from the alignments with the neighbour-join- 
ing method in MEGA 5.0 software, with Echinococcus multilocu- 
laris as the outgroup. Confidence intervals were obtained via 
1,000 bootstiap replications for each branch of the ttee. The 
identification of haplotypes and constmction of the haplotype 
network was undertaken by TGS 1.21 software using statistical 
parsimony [21]. The network estimation was mn at a 95% con- 
nection limit. Several sequences were obtained from GenBank 



Fig. 1 . Electrophoresis results of amplification product of cytb 
gene. Lane 1 , £ granulosus DNA; Lane 2, template free; Lane 3, 
marker; Lane 4, primers free. 

for similarity comparison and phylogenetic analysis. The com- 
plete mitochondrial genome of E. granulosus Gl (AF297617), 
G4 (AF346403), G5 (AB208063), G6 (AB235846) G7 
(AB23587), and G8 genotype (AB235848) were served as refer- 
ence sequences 122,23], and E. multilocularis (AB018440) as the 
outgroup [24]. Published mtiDNA sequences of cytb (AY278067 
and AB622276) [17,19] were used for comparison. 

PGR amplification products of cytb gene were analyzed by 
using 1% agarose gel electrophoresis, and the results showed 
tiiat a 1,609 bp band was amplified between 1,000-2,000 bp, 
which was consistent with the expected length of the target 
fragment (Fig. 1). DNA sequencing revealed that the cytb gene 
sequences of all 45 samples contained 1,068 nucleotides, and 
sequences were deposited in GenBank with accession numbers 
KG709522 to KG709566. The average contents of A, T, G, and 
G were 19.4%, 47.7%, 8.9%, and 24.0%, respectively revealing 
that AT content was greater than GC content. There were 11 
mutation sites among these sequences of cytb genes (Table 2). 
Only nucleotide substitutions were detected, while deletions or 
insertion mutations were not observed. High similarity was 
shown among all the sequences obtained in tiiis study (99.5% 
to 100%). In addition, all sequences in this study showed 
more than 88.0% identity with the previously published cytb 
sequences for E. granulosus genotypes. 

We deteaed 10 haplotypes within the 45 isolates obtained 
from the west of Ghina (Table 2). It was indicated that each of 
the haplotypes possessed intermediate host specificity. Haplo- 
type diversity (Hd) was 0.626, while nucleotide diversity (Pi) 
was 0.001. All the HI haplotype sequences shared 100% iden- 
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Table 2. Variation sites of cytb gene of Echinococcus granulosus from tine west of Cliina 



Haplotype 


Locality 


Host 


N" 


97 


120 


180 


285 


471 


680 


768 


850 


888 984 


1,038 


HI 


Qingliai, Sichuan, and Tibet 


Sheep Human 


27 


A 


A 


A 


T 


T 


0 


T 


C 


G C 


G 


H2 


Qinghai, Tibet 


Sheep 


5 




G 


















H3 


Qinghai 


Sheep 


3 


G 




G 










T 


T 




H4 


Qingtiai, Tibet 


Sheep 


4 






G 
















H5 


Qinghai 


Sheep 






















A 


H6 


Qinghai 


Sheep 












0 












H7 


Qinghai 


Sheep 










0 














H8 


Qinghai 


Sheep 




















A 




H9 


Sichuan 


Human 














T 










H10 


Tibet 


Sheep 
















C 









^Number of isolates. 
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Fig. 2. Neighbour-Joining phylogenetic tree of complete mito- 
cliondrial Echinococcus granulosus cytb sequences sinewing the 
phylogenetic relationship of haplotypes, H1 to H10. 




# Tibet O Sichuan ® Qinghai 

Fig. 3. Haplotype network for the 10 haplotypes (H1 to H10) 
identified among 45 isolates of Echinococcus granulosus collect- 
ed from western China. 



tity with a partial cytb gene sequence of a yak E. granulosus strain 
from China (AY278067) [17], and differed by only 1 nucleo- 
tide (a C/T substitution at the 263th site) from a cat E. gtmulo- 
sus strain found in lapan (AB622276) [19]. 

The phylogenetic tree showed that all isolates fr om the west 
of China were located on a single branch of the phylogenetic 
tree corresponding to E. granulosus genotype CI, and all haplo- 
types identified in this study were distinguished from each oth- 
er (Fig. 2). The haplotype network of 45 isolates of E. granulosus 
revealed that HI was the predominant haplotype (Fig. 3). HI 
emitted 7 additional haplotypes and a branch containing 2 
haplotypes (H3 and H4). 

To date, a few studies have explored in detail genetic diversity 
of E. granulosus based on cytb gene. The partial nucleotide se- 
quences of coxl (795 bp) and qftb (568 bp) have amplified 
from a presumptive echinococcal liver lesion of a yak in Sich- 



uan Province, and the results demonstrated that the yak was in- 
fected with E. granulosus CI [17]. Based on partial sequences of 
cytb gene (549 bp) and rmL (570 bp), Li et al. [18] reported that 
33 of 53 isolates obtained from humans in Sichuan and Qing- 
hai were E. gmnulosus CI, while the remaining 20 isolates were 
E. multilocularis. In addition, cytb (1,068 bp) and coxl (1,069 
bp) gene sequences revealed that the causative agent of CE in a 
cat from Russia was E. granulosus s.s. (CI) [19]. In our study, us- 
ing complete sequences of cytb gene (1,068 bp) we identified 
all 45 isolates obtained from westem China as belonging to the 
strain of E. gmnulosus genotype CI. The data obtained in tiiis 
study are more representative and comprehensive than the par- 
tial sequence. The complete sequence may provide a more ef- 
fective evidence to verify the validity of the genotypes of E. 
granulosus. 

In recent years, studies based on nuclear and mitochondrial 
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genes indicated that the Gl genotype is the predominant epi- 
demic strain in the west of China, while the G3 and G6 strains 
also exist. Based on coxl (789 bp) and efla (656 bp), the iso- 
lates of E. granulosus s.s. (n = 181) and E. canadensis G6 (n = 1) 
were identified successfially in Qinghai, Sichuan, and Xingjiang 
Province [2]. Phylogenetic analysis based on nadl (332 bp) 
and coxl (450 bp) revealed that 28 isolates of E. granulosus be- 
longed to G1-G3 complex (£. granulosus s.s.), and 2 isolates 
were placed in G6-G10 complex [E. canadensis) in Qinghai [11]. 
Yan et al. [13] reported that 84 isolates of E. granulosus from the 
Tibetan plateau in China were sequenced for the complete 
nadl (894 bp) and atp6 (513 bp) gene, and the results showed 
that 82 isolates were identified as Gl genotype, while 2 as G3 
genotype. Our findings are similar to previous results that Gl 
genotype of E. gtmulosus is the main epidemic genotype in the 
west of China, but we did not detect G3 or G6 genotype in 45 
samples of E. granulosus. 

Our results showed low genetic diversity among the 45 iso- 
lates of E. gmulosus fr-om the west of China. In addition, the 
similarities of cytb genes of 45 isolates was very high (99.5- 
100%), indicating that the genetic diversity among the 45 iso- 
lates was small. 

In conclusion, Gl genotype was identified to be the only 
genotype in all 45 samples of E. gmnulosus from Tibet, Sich- 
uan, and Qinghai provinces. The complete cytb gene investi- 
gated in this study enriched the genetic information among E. 
granulosus genotypes from the west of China. 
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